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ABSTRACT 

We present Spitzer infrared imaging of the peculiar galaxy pair Arp 107, and compare with 
an optical Ha map and a numerical model of the interaction. The [3.6] — [4.5] colors of clumps 
in the galaxy do not vary around the ring-like primary spiral arm and are consistent with those 
of stars, thus these bands are dominated by starlight. In contrast, the [5.8 /imj — [8.0 ^m] colors 
are consistent with those of interstellar dust, and vary by about 0.2 magnitudes around the 
ring/spiral, with redder colors associated with regions with stronger star formation as indicated 
by Ha and mid-infrared luminosity. The [4.5 /im] — [5.8 ^m] colors for clumps in this arm are 
bluer than dust and redder than stars, and vary by 1.3 magnitudes around the arm. This color 
is therefore a measure of the relative number of young stars to old stars, with a redder color 
indicating a higher proportion of young stars. There is a clear azimuthal sequence in the [4.5] — 
[5.8] color around the arm, indicating a sequence in average stellar age. The L/fa/Lg o^im ratio 
varies around the arm by a factor of ~7; this variation may be due to extinction or to PAH 
excitation by non-ionizing photons. Our model pf Arp 107 accounts for the general morphology 
of the system, and explains the age variation along the arm as the result of differences in the time 
of maximum compression in the arm. Using Spitzer colors, we are able to distinguish background 
quasars and foreground stars from star forming regions associated with Arp 107. 

Sii.hip.r.f. hp/],din.n.i: ecalaxifis: starbnrsts — eralaxifis: interar.tinns — eralaxifis: inHiviHnal (Am 1fl7V 



1. Introduction 

Tidal interactions, collisions, and mergers of 
galaxies play a major role in the evolution of galax- 
ies, drastically modifying both the morphology 
and the star formation rates of galaxies (see re- 
view by Struck 1999). Interaction-induced star 
formation was first investigated by Larson & Tins- 
ley (1976) and Struck-Marcell & Tinsley (1978), 
who provided evidence that interacting galaxies 
have a larger scatter in star formation proper- 
ties than isolated galaxies. The Infrared Astro- 
nomical Satellite (IRAS) revealed a new popula- 
tion of high infrared luminosity galaxies with large 
star formation rates (Soifer et al. 1987; Smith et 
al. 1987) that are the result of mergers (Sanders 
et al. 1988). Mid- infrared observations with the 
Infrared Space Observatory (ISO) suggested that 
the mid-infrared is a good measure of the global 
star formation rate in galaxies (Roussel et al. 
2001; Forster Schreiber et al. 2004). The mid- 
infrared emission is dominated by broad dust emis- 
sion features, which are believed to be due to 
polycyclic aromatic hydrocarbons (PAHS) (Leger 
& Puget 1984; AUamandola, Tielens, & Barker 
1985). With the advent of the Spitzer infrared 
telescope (Werner et al. 2004), higher angular res- 
olution mid-infrared imaging of galaxies is now 
possible, making feasible the detailed study of in- 
dividual star formation complexes in interacting 
systems. 

In this paper, we present Spitzer infrared im- 
ages of the interacting galaxy pair Arp 107 (UGC 
5984), and compare with an optical Hct map to 
investigate the properties of star forming clumps 
with the galaxy. We also present a new numerical 
model of the interaction. The larger galaxy in the 
pair has a Seyfert 2 nucleus (Keel et al. 1985), a 
strong primary spiral arm with a partial ring-like 
structure, and a tidal tail (Arp 1966; sec Figure 
1). It is connected to the smaller companion by a 
bridge. Arp 107 is at a distance of 138 Mpc (Keel 
et al. 1985), asssuming Ho = 75 km Mpc"-'^. 

Globally, Arp 107 has a relatively low nor- 
malized star formation rate compared to other 
interacting galaxies. Its far-infrared luminosity 
(40 - 120 nm) is 1.1 x 10^° Lq, similar to 

^Astronomy Department, Cornell University, Ithaca NY 
14853 and Chercheur Associe, Observatoire de Paris, F- 
75014, Paris, Prance 



the median for other Arp Atlas galaxies, how- 
ever, it has a far- infrared to blue luminosity ra- 
tio Lfir/I^b = 1-6 and an IRAS 60-to-lOO /zm 
flux ratio Fgopm/^'ioo/im — 0.19, at the low end 
of the range for Arp galaxies (Bushouse, Lamb, & 
Werner 1988). Arp 107 is relatively rich in gas, 
with a total HI mass of 1.2 x 10^° Mq (Bushouse 
1987) and a molecular gas mass of 9.6 x 10^ Mq 
in a 22" beam (Zhu et al. 1999), assuming the 
standard Galactic Ico/^Hi ratio (Bloemen et al. 
1986). 

2. Observations 

2.1. Spitzer Observations 

Images of Arp 107 were made on December 
17, 2004 in the 3.6, 4.5, 5.8, and 8.0 nm broad- 
band filters of the Spitzer Infrared Array Camera 
(IRAC; Fazio et al. 2004). The full bandwidths of 
these filters at 50% of the average in-band trans- 
mission are 0.75 /um, 1.01 iim, 1.42 /Lim, and 2.93 
Hm, respectively (Fazio et al. 2004). A total of 23 
dithered exposures of 12 seconds each were made 
per filter, giving a total integration time per filter 
of 276 seconds. The standard post-pipeline Ba- 
sic Calibrated Data (post-BCD) mosaicked images 
were used for this analysis, providing a platescale 
of 1'.'2 per pixel. The FWHM of the point spread 
function at 3.6, 4.5, 5.8, and 8.0 fim are 1'.'7, 1'.'7, 
1'.'9, and 2'.'0, respectively (Fazio et al. 2004). The 
total field of view mapped with IRAC is 9'.0 x 
14.'8, with the long axis at a position angle of 19.8° 
east of north. For 3.6, 4.5, 5.8, and 8.0 /xm, the rms 
noise levels in the final maps are 0.0072, 0.0099, 
0.041, and 0.041 MJy sr"\ respectively. 

Additional observations of Arp 107 were made 
on December 3, 2004, using the 24 fim filter (band- 
width « 5 /im) of the Spitzer Multiband Imaging 
Photometer (MIPS; Rieke et al. 2004). Two 10- 
second cycles were made with the standard pho- 
tometry AOT, giving an effective exposure of 312 
seconds. The pixel size of the post-BCD images is 
2'.'45, and the FWHM of the point spread function 
is 6". The total field of view mapped is 7.'4 x 8.'1, 
with the long axis of the field of view at a position 
angle of 23.9° east of north. The final rms noise 
level in this map is 0.048 MJy sr~^. 



2 



2.2. Optical Observations 

Optical images of Arp 107 were obtained with 
the 0.6 meter Erwin Fick telescope in Boone Iowa 
on 2005 April 8, 14, and 15, under clear dark skies. 
An 1024 X 1024 Andor DU434-BV CCD was used 
with a pixel size of I'.'OS, giving a field of view 
of 17'.9. A total of seven 120 second exposures 
were made in a broadband R filter, along with 7 
X 450 second images in a redshifted Ha filter cen- 
tered at 6814A with FWHM 73A. For Arp 107, 
this filter contains both Ha and the [N 11] A6583 
line. The data were reduced in the standard way 
using the Image Reduction and Analysis Facility 
(IRAF^) software. Continuum subtraction was ac- 
complished using the scaled R band image. The 
standard star Feige 56 was observed for Ha cal- 
ibration. The Ha calibration is estimated to be 
accurate to ~30%. 

3. RESULTS 

3.1. The Optical Images 

The final co-added Fick R band image of Arp 
107 is presented in Figure 1. In this image, the 
strong primary arm/ring- like structure is appar- 
ent in the south, which connects to a tidal tail to 
the northwest. A bridge connects the main galaxy 
to a compact companion to the northeast. The 
bridge appears double in this image, and a plume 
is visible to the northeast of the companion. 

The continuum-subtracted Ha-|-[N II] map is 
shown in Figure 2. Seven star formation regions 
were detected in the ring/arm, in addition to the 
Seyfert nucleus. The Ha-|-[N II] luminosities of 
these HII region complexes range from »4 x 10^^ 
erg s~^ to 3 x 10^° erg in 5" circular aper- 
tures (see Table 1). These are similar to those 
of the most luminous HII regions in ring galaxies 
(Marston & Appleton 1995; J. Higdon 2004, pri- 
vate communication). The total Ha-|-[N II] lumi- 
nosity for this system, excluding the Seyfert nu- 
cleus, is «1.1 X lO^'^ erg s"-*^, at the lower end 
of the range for the ring galaxies in the Marston 
& Appleton (1995) sample. For the Seyfert nu- 
cleus, ljHa+[N II] ^ 2 X 10*^ erg s~^, consistent 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which arc operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



with the Keel et al. (1985) measin-ements. The 3(t 
'^Ha+{N II] upper limit for the companion is w 4 x 
10^° erg s~^. Approximately correcting for [N II] 
in the filter and excluding the Seyfert nucleus, the 
far-infrared to Ha luminosity ratio for Arp 107 is 
^fir/^Hoi ~ 520. This is typical of interacting 
galaxies (Bushouse 1987). 

3.2. The Spitzer Images 

The Spitzer images at the five wavelengths are 
presented in Figure 3. The morphology of this 
system varies with wavelength. The bridge is clear 
at 3.6 (Um but almost disappears by 8.0 fjbm. In 
contrast, the knots in the northwestern part of the 
ring/arm are more prominent at 8.0 /xm than at 
3.6 jum. The companion is very faint at 24 fjbm but 
the ring/arm is clearly seen. The double bridge 
structure is visible in the two shorter wavelength 
IRAC bands. 

At 3.6, 4.5, 5.8, 8.0, and 24 /xm, the total flux 
densities for the larger galaxy in the pair are 21.7 
mJy, 14.2 mJy, 19.9 mJy, 38.0 mJy, and 53.4 mJy, 
respectively. For the smaller galaxy, the total flux 
densities are 13.2 mJy, 8.5 mJy, 13.0 mJy, 4.4 
mJy, and 0.48 mJy, respectively. Absolute calibra- 
tion uncertainties are <10% (IRAC Data Manual; 
MIPS Calibration Manual). Statistical uncertain- 
ties are <1%, except for the companion at 8.0 and 
24 /xm (2% and 17%, respectively). Ambiguities 
in selecting the boundaries of the galaxies and sky 
subtraction introduce an additional uncertainty of 
<20% in most cases, 40% at 5.8 ^,m, and 60% for 
the companion at 24 /im. 

The global Lna/Ls.o/^m and hHa/'^^..&^lm ra- 
tios are «0.0083 and 0.012, respectively, for the 
main galaxy (excluding the Seyfert nucleus), 
and <0.025 and <0.0061 for the companion. 
For comparison, the global log l^Ha I^^.Q^im and 
Ljfa/Ls.e^im ratios for the spiral galaxy NGC 7331 
are f»0.013 and 0.030, respectively, calculated us- 
ing the Young et al. (1996) Ha measurement and 
the Spitzer flux densities of Regan et al. (2004). 

3.3. IRAC Colors of Star Forming Clumps 

In the 8 ^m image, 29 clumps were identified by 
eye in the vicinity of Arp 107. These are labeled 
on the 8 /xm image in Figure 4, encircled by 5" 
radii regions, and are listed in Table 1. Note that 
some fainter clumps are not included. When a 
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second fainter clump was very close to a brighter 
clump so that their regions strongly overlapped, 
the second clump was excluded. 

Photometry of these clumps was done using 
the IRAF daophot routine, with 4.1 pixel (5") 
radii apertures (Table 1). For background sub- 
traction, the local galaxian background was deter- 
mined using the mode in an annulus surrounding 
the source, with an inner radius of 5 pixels (6") 
and an outer radius of 10 pixels (12"). To estimate 
the uncertainty in the Spitzer colors of the clumps 
due to background subtraction, we experimented 
with varying the method of determining the sky. 
In addition to the statistical uncertainties deter- 
mined from the rms in the background annuli, in 
calculating the colors we added in quadrature a 
second uncertainty term, determined from com- 
paring the clump colors obtained with the above 
method with those obtained with slightly larger 
inner and outer annulus radii of 7 pixels (8'.'4) and 
12 pixels (14'.'4). This additional term increases 
the median size of the final uncertainties in the 
[3.6] - [4.5], [4.5] - [5.8], [5.8] - [8.0], and [3.6] 
- [8.0] colors by 15%, 33%, 8%, and 66%, respec- 
tively. 

No color corrections were made in determining 
the IRAC colors. Aperture corrections were in- 
cluded, although they have only a tiny effect on 
the IRAC colors of these clumps. As given in the 
IRAC Data Manual, the aperture corrections are 
1.06, 1.06, 1.07, and 1.09 for IRAC bands 1-4, 
thus the changes to the IRAC colors due to the 
aperture corrections are negligible. These clumps 
are unresolved or barely resolved at the resolu- 
tion of the mosaicked images (2'.'4 — 3'.'6, com- 
pared with 2'.'0 — 2'.'4 for bright point sources in 
the field). 

In Figures 5 and 6, the IRAC colors of the 
brightest 8 jiui clumps in Arp 107 are plotted 
(black crosses) and identified. Both statistical and 
background selection uncertainties arc included in 
these plots, but absolute calibration uncertainties 
(<10%; IRAC Data Manual) are not included. 
Figures 5 and 6 also contain the colors of MOIII 
stars (open dark blue square) (from M. Cohen 
2005, private communication) and the mean colors 
for the field stars in Whitney et al. (2004) (ma- 
genta triangle). The colors of normal stars all 
lie within «0.5 magnitudes of 0,0 on these plots 
(M. Cohen 2005, private communication). The 



green diamonds in these figures mark the predicted 
IRAC colors for interstellar dust from Li & Draine 
(2001) (see Table 5 in that paper), for a very large 
range of interstellar radiation field (ISRF) intensi- 
ties (from 0.3 — 10,000 x that in the solar neigh- 
borhood). Note that the IRAC colors of the Li & 
Draine (2001) dust model varies very little with 
ISRF. The red circles in Figure 5 show the loca- 
tions of the Sloan Digitized Sky Survey quasars in 
the Spitzer Wide- Area Infrared Extragalactic Sur- 
vey (SWIRE) Elais Nl field (Hatziminaoglou et al. 
2005). These quasars have redshifts between 0.5 
and 3.65; since their spectral energy distributions 
are power laws, their infrared colors do not vary 
much with redshift. 

The colors of Source #23 (the bright source to 
the west of the galaxy) are similar to those of the 
quasars. Thus this is likely a background quasar, 
unassociated with Arp 107. The nucleus of the 
companion, source #28, is very blue, with stel- 
lar colors at all wavelengths. Thus it contains al- 
most all old stars and little ISM. This is consis- 
tent with its absorption-line optical spectrum (W. 
Keel 2005, private communication). The nucleus 
of the primary galaxy, #14, is also quite blue, sug- 
gesting an old stellar population, but not as old as 
the companion nucleus #28. This nucleus also has 
some interstellar matter, as shown by its large [5.8] 
— [8.0] color. Source #9, a concentration inside 
the radius of the ring, also has colors indicative 
of an old star, suggesting it is a foreground star. 
This is confirmed by its optical spectrum, which 
is that of an M1-2V star (W. Keel 2005, private 
communication) . 

The clumps in the upper part of Figure 6 have 
[5.8] — [8.0] colors close to that predicted by the 
Li & Draine (2001) dust model, indicating that 
these wavcilcmgths arc; dominated by dust. The 
8.0 /im IRAC band includes the strong 7.7 and 8.6 
yLtm PAH emission features as well as a dust con- 
tinuum from transiently-heated very small grains, 
while the 5.8 /im IRAC band contains the 5.7 fiva 
PAH feature (Li & Draine 2001). In contrast, the 
[4.5] — [5.8] colors of the Arp 107 clumps are much 
bluer than those of dust, thus there are stellar con- 
tributions at 4.5 /xm. The [3.6] — [4.5] colors of 
these clumps are similar to those of stars (Figure 
5), indicating that these bands are dominated by 
starlight. 

In Figures 7 and 8, we have plotted the IRAC 
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colors of the fainter Arp 107 clumps. The colors 
of Source #19, off the western side of the galaxy, 
indicate that it is likely a background quasar. 
Sources #1, #2, and #3 (south of the ring), #12 
(west of the ring), #27 (in the bridge), #29 (at 
the northern tip of the tail), and #25 (northwest 
of the ring) have [3.6] - [4.5] and [4.5] - [5.8] col- 
ors like old stars, but have excess emission in the 
8.0 /xm band compared to old stars. This suggests 
that either these are very old star formation re- 
gions associated with the galaxy, with some PAH 
emission but mainly cooler stars, or they are fore- 
ground stars, and the 8 /im band is contaminated 
by extended dust emission from Arp 107. 

3.4. An Azimuthal Color Gradient Around 
the Ring/ Arm 

In the [4.5] - [5.8] vs. [5.8] - [8.0] plot for the 
bright clumps (Figure 6), there is a sequence of 
sources at an almost constant [5.8] — [8.0] color of 
ps2, with increasing [4.5] — [5.8] colors. The se- 
quence from left to right at the top of Figure 6 is 
a sequence in increasing ratio of young stars/old 
stars. This sequence, #4, #5, #7, #10, #21, 
and #26, is the order these clumps appear in 
the ring/ arm, going in a counterclockwise direc- 
tion from #4 to #26. The northern part of the 
ring/arm (clumps #21 and #26) has a younger 
stellar population than the southern part (#4 and 
#5) on average. This trend in the [4.5] — [5.8] 
color is better shown in Figure 9, plots of position 
angle (east of north) vs. the IRAC colors. There 
is a strong azimuthal variation in the [4.5] — [5.8] 
color of pal. 3 magnitudes. In contrast, no vari- 
ation is visible in the [3.6] — [4.5] color around 
the ring/arm. The approximately constant [3.6] 
— [4.5] color of Ri 0.0 around the ring/arm sug- 
gests that these bands are dominated by emission 
by cool stars. 

In the [5.8] — [8.0] color, the variation is small 
(0.2 magnitudes) compared to that in [4.5] — [5.8], 
and the northern clumps are redder in [5.8] — [8.0] 
than the southern clumps. This small variation 
and the similarity to the colors of interstellar dust 
indicates that these two bands are dominated by 
dust emission, and the dust spectrum does not 
vary strongly around the ring/arm. The [5.8] — 
[8.0] color increases in the north, in the regions 
where the Ha and mid-infrared luminosities get 
larger (see Section 3.6). Thus this color gets red- 



der as the ISRF gets stronger, as predicted by the- 
oretical models of interstellar dust (see Table 5 in 
Li & Draine 2001). As the ISRF increases, the un- 
derlying mid-infrared continuum from very small 
grains increases, contributing to the 8 /xm IRAC 
band, and producing redder [5.8] — [8.0] colors. 

Figure 9 shows a large variation in the [3.6] — 
[8.0] color around the ring/arm (bottom panel), 
which is mainly due to changes in the [4.5] — [5.8] 
color (second panel). 

3.5. The IRAC-MIPS Colors of Clumps 

Prom the MIPS 24 /xm map, 24 /xm fluxes were 
measured using a 5" radius aperture. The back- 
ground was determined using the mode in an an- 
nulus with an inner radius of 10" and an outer ra- 
dius of 15". The 5" radius aperture is too small to 
include all of the light from the clump, given the 
w6" FWHM PSF for MIPS 24 /zm observations 
and the 2.45"/pixcl plate scale. Using a larger 
aperture radius, however, is not a viable option, 
since the clumps are crowded. Thus a substantial 
aperture correction is required. According to the 
MIPS Data Handbook, a 5" aperture includes only 
~50% of the total light from a point source. Test- 
ing this using large aperture data for the isolated 
source #23 (which may be a background quasar; 
see Section 3.3) and two bright sources in the field 
far from the galaxy, we get consistent aperture 
corrections. This correction and the 8 fiui aper- 
ture correction given above increases the [8.0] — 
[24] colors by 0.65 magnitudes. No color correc- 
tions were included in determining these colors. 
To estimate the uncertainty in the colors due to 
the selection of background regions, we also cal- 
culated colors with a smaller background annulus 
with inner radius of 5" and outer radius of 10", 
and compared with the original colors. This addi- 
tional color uncertainty was added in quadrature 
to the statistical uncertainty. This extra factor 
increases the sizes of the errorbars on the IRAC- 
MIPS colors of the clumps by a factor of f«2.7 over 
the statistical errors. 

In Figure 10, the [5.8] - [8.0] vs. [8.0] - [24] 
aperture-corrected colors for the clumps are plot- 
ted. The green diamonds show the locations of 
the Li & Draine (2001) dust models for various 
ISRFs. The red circles are the locations of the 
Hatziminaoglou et al. (2005) quasars. Stellar pho- 
tospheres are expected to lie near [5.8] — [8.0] » 
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0, [8.0] - [24] « 0. 

The foreground star #9 again has colors simi- 
lar to those of stars, as expected. The [8.0 /um] — 
[24 /Ltm] color of Source #23, like its IRAC colors, 
is similar to those of the quasars, further support- 
ing the hypothesis that it is a background quasar. 
Sources #1 and #12 are also in the vicinity of the 
quasars. In this plot, the Seyfert nucleus, clump 
#14, has colors similar to those of quasars, with 
a very red [8.0] — [24] color. The nucleus of the 
companion, #28, shows a small excess relative to 
the expected [8.0] — [24] color of stars («0.0), thus 
there may be a small amount of cold dust associ- 
ated with this nucleus. 

The [8.0] — [24] colors of the brightest ring/arm 
clumps, #10, #21, and #26, are similar to those 
expected for standard dust models with weak ra- 
diation fields (approximately solar neighborhood). 
This is consistent with their relatively low ob- 
served average 5 — 8.5 surface brightnesses 
of PslO Lq pc~^, which corresponds to an ISRF 
of «2 times that in the solar neighborhood, ac- 
cording to Forster Schreiber et al. (2004). In Fig- 
ure 11, the [8.0] — [24] colors for the ring/arm 
clumps are plotted as a function of position an- 
gle. A possible trend is seen, in that the clumps 
near the southern part of the ring/arm, #4 — #7, 
have slightly bluer [8.0] — [24] colors than those 
in the northwest, however, this is uncertain. If 
this is confirmed, it would imply that the ISRF is 
stronger in the northern part of the ring/arm, as 
also indicated by the IRAC colors. 

3.6. Ha to Mid-Infrared Ratios 

Of the extra-nuclear star formation regions in 
Arp 107, only clumps #4, #5, #7, #10, #16, 
#21, and #26 were detected in Ha, with clumps 
#4 and #5 detected at the 5<j level (Section 3.1; 
Table 1). In the top panel of Figure 12, for the 
selected clumps we compare the ratio of the Ha 
hmiinosity L/ja to the 3.6 /xm luminosity La.g/^m 
with the 8.0 fj,m magnitude, correcting for [N II] 
in the Ha filter assuming [N II] /Ha « 0.3. In the 
second panel of Figure 12, we plot L^/a/Ls-O/im vs. 
the 8.0 iim magnitude. For the clumps detected 
in Ha, log Li/a/Ls.o/im « -1-3 to -2.2. 

Figure 12 shows that clumps #4, #5, #6, 
and #7 arc deficient in observed Ha compared 
to clumps #10, #21, and #26. The ratios 



Lffa/Ls.OMm and Lna/Ls.e^m vary systemati- 
cally around the ring/arm (Figure 13), with 
Lffa/Ls.OAtm ~ 2 — 11 times larger in the north 
and west compared to the south and east. 

One possible explanation for the low Lija/Lg.o/iim 
ratios in clumps #4 — 7 is that these clumps have 
some contributions to their 8 fim emission from 
starlight. However, the fact that their [5.8] — 
[8.0] colors are similar to those of interstellar mat- 
ter and also to the other clumps in the ring/arm 
indicates that this is not the case. 

A second hypothesis is that in clumps #10, 
#21, and #26, strong radiation fields may be 
destroying the PAH molecules contributing to 
the 8.0 /xm Spitzer data, causing an increase in 
Ljja/Lg.Ojum- Using ISO data for a large sample 
of spiral and starburst galaxies, Forster Schreiber 
et al. (2004) found a good correlation between 
the mid-infrared luminosity and inferred star for- 
mation rates for 5 — 8.5 /im surface brightnesses 
below RilO^ Lq pc~^; above this surface bright- 
ness, the mid-infrared appeared weaker than ex- 
pected, suggesting PAH destruction. In Arp 107, 
however, such large-scale PAH destruction is not 
expected to happen, since the radiation fields are 
not intense enough. For the brightest clumps in 
the Arp 107 ring/arm, the 5 — 8.5 /xm surface 
brightnesses of «10 Lq pc~^ are well below this 
threshold, thus large-scale PAH destruction is not 
expected. Furthermore, the [5.8] — [8.0] colors 
of clumps #21 and #26 are slightly redder than 
those of #4 and #5, opposite to what would be 
expected if the 8 /xm broadband flux is being de- 
pressed due to destruction of PAHs by an intense 
UV field. 

A third possibility is that non-ionizing photons 
in the southern part of the ring/arm contribute 
significantly to heating the PAH molecules, and 
therefore the 8 /xm band is not a perfect tracer 
of star formation. PAH excitation by non-ionizing 
photons has been noted in reflection nebulae (Sell- 
gren et al. 1990; Uchida, Sellgren, & Werner 1998; 
Uchida et al. 2000; Li & Draine 2002), and has 
been suggested contributor to the global 

mid-infrared of spiral galaxies based on ISO data 
(Boselh, Lequeux, & Gavazzi 2004). The PAH fea- 
tures may therefore be a tracer of B stars as well 
as O stars (Peeters, Spoon, & Tielens 2004). If 
the southern and eastern part of the ring/arm con- 
tains a significant post-starburst population, these 
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modcratcly-young stars may be sufficient to excite 
the PAH molecules and produce the observed mid- 
infrared emission. The bluer [4.5] — [5.8], [5.8] — 
[8.0], and [8.0] — [24] colors of the southern clumps 
#4 and #5 compared to clumps #21 and #26 sug- 
gest that the mid-infrared continuum from very 
small grains is weaker in the south, and therefore 
the exciting UV field is weaker. This is consis- 
tent with the hypothesis of an older stellar popu- 
lation in the south. This possibility should be in- 
vestigated further with optical and infrared spec- 
troscopy and stellar population synthesis model- 
ing. 

A fourth possibility is that the Ha emission 
is more extincted in the clumps in the southern 
and eastern part of the ring/arm compared to the 
northwest. A relative extinction of Ay « 1.5^0 g 
magnitudes around the ring/ arm could account for 
the observed deficiency in the Ha luminosity. This 
possibility should be tested further with optical 
spectroscopic measurements of the Balmer decre- 
ment and/or high angular resolution CO and HI 
maps. This will determine whether the observed 
azimuthal trend in IjHa/^&.Qfim in the ring/arm 
is due to extinction or to PAH excitation by non- 
ionizing photons. 

In the L^fa/La g/^m vs. [8.0] plot, the same 
trends are seen: the 'LHal^3.& ratios for clumps 
#4, #5, #6, and #7 are lower than clumps #10, 
#21, and #26. The foreground star #9 has a very 
low Lij„/L3.6jum ratio, as expected. 

4. A Numerical Model of the Encounter 

The morphology of the primary galaxy in Arp 
107 consists of a small central bar-like bulge and 
a single strong spiral arm with a ring-like shape. 
Most of this arm is quite circular, though the circle 
is not centered on the bulge and nucleus. Given 
this, and the fact that there is material between 
the companion and the primary, like parts of a 
bridge, the system also looks like an asymmet- 
ric colliding ring galaxy (Lynds & Toomre 1976; 
Theys & Spiegel 1977; Appleton & Struck 1996). 
On the other hand, off-center ring waves usually 
look like ovals, not one-armed spirals (Toomre 
1978; Struck-MarceU 1990; Smith & Wallin 1992). 
In addition, our Spitzer maps reveal a strong az- 
imuthal gradient in the ratio of young stars to old 
stars around the ring of Arp 107. This asymme- 



try is also consistent with the hypothesis that this 
structure was produced by an off-center collision 
between two galaxies. Theoretical models suggest 
that such encounters lead to strong azimuthal vari- 
ations in the gas density and star formation rate in 
collisional rings (Appleton & Struck-MarceU 1987; 
Struck-MarceU & Appleton 1987; Charmandaris, 
Appleton, & Marston 1993). 

No complete velocity field of Arp 107 is cur- 
rently available. However, Keel (1996) obtained a 
long-slit optical spectrum of Arp 107 along a posi- 
tion angle of 41° east of north. He found that the 
southwestern side of the main disk is redshifted 
compared to the nucleus. The companion has a 
velocity cz = 10,572 ± 50 km s~^, redshifted by 
200 km s~^ from the main galaxy (W. Keel 2005, 
private communication). 

To resolve the nature of this collisional system 
and to better understand the structure of the star 
formation in the arm we have attempted to numer- 
ically model it. Given that we have limited kine- 
matical information and very little information on 
the gas distribution, we have not attempted to 
model it in great detail, but rather to produce a 
reasonable preliminary model. 

4.1. Details of the Model 

We have used the SPH code described in Struck 
(1997). The version used here is essentially the 
same as that used in Struck & Smith (2003). 
Briefly, this code uses rigid dark halo potentials 
for the two galaxies, and hydrodynamic forces are 
computed on a grid with fixed spacing. Local 
gravitational forces are computed between parti- 
cles in adjacent cells. The galaxy disks are made 
up of gas particles and coUisionless star particles of 
equal mass. In these models the particle numbers 
were as follows: 34,550 primary disk gas particles, 
18,090 companion disk gas particles, 7890 primary 
star particles and 2490 companion star particles. 

Because the companion appears to be an early 
type galaxy, its star particles were initialized with 
a large range of distances from the central plane 
and a large velocity dispersion perpendicular to 
that plane to give a very thick disk or bulge-like 
structure. A substantial gas disk was also added to 
the companion, even though there is no evidence 
for the existence of a gas disk. The goal was to 
see how it would be affected by the collision, and 
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whether we would expect it to be detectable in 
current or future observations. 

We adopt the following scaling constants for the 
model: time unit = 400 Myr, length unit = 2.0 
kpc, and mass unit = 1.3 x 10^° Mq. With these 
scalings we obtain a peak rotation velocity in the 
primary of about 170 km/s, which is reasonable 
within the very limited constraints. 

The mass ratio of the two galaxies, determined 
by comparing the mass of each within a radius of 
10 kpc, is 0.16. The form of the rigid potential 
used gives a test particle acceleration of, 



where is a halo mass scale, e is a core radius 
(set to 4.0 and 2.0 kpc for the primary and com- 
panion, respectively), and the index specifies 
the compactness of the halo. For the primary we 
use n/, = 1, which gives a flat rotation curve at 
large radii. For the companion we take = 1.35, 
which gives a declining rotation curve at large 
radii. This latter choice seems to produce some- 
what better results than two flat rotation curves, 
but we have not undertaken a detailed examina- 
tion of the effects of different values for these ex- 
ponents. Because the two galaxies have different 
halo potentials, the effective mass ratio (i.e., the 
ratio of masses contained within radii equal to the 
galaxy separation) is large for small separations, 
and small for large large separations. This effect 
makes encounters between the two galaxies quite 
impulsive. The model presented here does not in- 
clude the effects of dynamical friction, but other 
runs included a Chandrasekhar-like frictional term 
(see Struck & Smith 2003). This term was found 
to be negligible except for brief times near closest 
approach. It had a small effect at the first closest 
approach, but ultimately draws the two galaxies 
into a merger at times past the present (the sec- 
ond close encounter). 

The initial disk sizes were as follows: primary 
gas disk =15 kpc, primary star disk =12 kpc, 
companion gas disk = 9.0 kpc, and companion star 
disk = 3.6 kpc. These are relatively arbitrary val- 
ues until 21 cm mapping of the HI distribution 
becomes available. The companion disk was ini- 
tialized in the computational x-y plane. The pri- 
mary disk was initialized in the x-y and then tilted 



20° around the x-axis such that positive y values 
had positive z values. After the run the viewing 
angle was adjusted by rotating the system —50° 
around the x-axis, 30° around the y-axis, and 90° 
around the z-axis. The initial position and velocity 
vectors of the companion relative to the primary 
center were: (7.0 kpc. 7.0 kpc, 8.0 kpc) and (70 
km s~^, -44 km s~^, -24 km s~^). 

4.2. Model Results 

The dashed lines in the top panels of Figure 
14 show that, although the companion collides 
with the outermost part of the primary gas disk, 
it then flics a considerable distance radially out 
from the disk center and below it (assuming it be- 
gan 'above' it, i.e., at positive z coordinate values). 
The model seems to fit the observed morphology 
best at a time just before or after the companion 
again passes through the primary disk. This is 
in accord with the significant line-of-sight velocity 
difference between the two galaxies. 

The lower panels in the figure show that at ear- 
lier times there was a siibstantial bridge between 
the two galaxies, and the primary had much more 
of the two-armed appearance characteristic of pro- 
grade encounters. However, by the present time, 
most of the bridge material has fallen onto one or 
the other galaxies. What little is left is projected 
along the line-of-sight, and not very visible. Some 
bridge gas originating from the companion con- 
tinues to fall onto the primary. This mass transfer 
was greater in the past, and could plausibly play a 
role in feeding the Seyfert nucleus of the primary. 

We have run a number of additional exploratory 
models, and in many of them the companion has 
a substantial tidal tail. Depending on initial and 
orbital parameters these tails can have many dif- 
ferent orientations. In some cases tail material 
might account for the debris to the east of the com- 
panion and north of the primary. In the present 
model, however, debris to the east and west of the 
companion consists mostly of stars pulled off the 
primary disk (see Figure 14). 

In these models the primary arm originates as 

a tidal counter-tail drawn from primary disk ma- 
terial on the opposite side from the companion's 
closest approach point. The bridge arm originates 
on the companion side near the beginning of the 
models, and its outer part becomes the bridge as 
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described above. The inner part of the bridge arm 
is overrun by the primary arm at about the present 
time in many of the models. The models suggest 
that the pccuUar 'nob' on the eastern side of tlie 
primary arm is a remnant of the bridge arm. In 
sum, the models suggest that the one-armed ap- 
pearance of the primary galaxy is the result of the 
asymmetric evolution of an essentially two-armed 
disturbance. This asymmetric evolution is the re- 
sult of substantial m=0 (ring-like) and m=l (per- 
turbed center) components of that disturbance. 

The new Spitzer observations resolve many 
star-forming clumps in the primary arm. Inter- 
estingly, there seem to be clumps along almost 
the whole arm, and the star formation is not sig- 
nificantly stronger in the inner arm regions where 
one might expect the gas density to be greater. 
Moreover, as described above, there is evidence 
that clump populations are younger as one tracks 
the arm outward and around. 

The model results in the lower panels of Fig- 
ure 14 offer some possible explanations for these 
trends. These show that the arm has narrowed 
and stretched in the last 400 Myr. The overlaps 
of test particle rings show that compressions in- 
creased greatly on the leading edge of the arm 
in the west in the first half of this time period. 
The figure shows that compressions in the outer- 
most parts of the arm develop later. The feedback 
algorithm in the models shows that star-forming 
clumps appear rather stochastically at late times 
in the outer tail, which is consistent with the test 
particle results convolved with a random element 
in the way gas densities push over a threshold. We 
suggest that this may be how chimps #16, 21 and 
29 formed. Alternatively, clumps #21, 22, 24 and 
26 might be the result of companion debris inter- 
acting with the primary disk. 

We should emphasize that this relatively straight- 
forward picture of the nature of the star formation 
in the primary arm is strongly based on the Spitzer 
infrared colors. The situation would be consider- 
ably murkier if we had only Ha data to compare 
to the models. 

5. Conclusions 

We see a strong azimuthal variation of 1.3 mag- 
nitudes in the [4.5] — [5.8] colors of clumps in the 
ring/primary arm of Arp 107. These colors are in 



between those of stars and dust, implying an az- 
imuthal sequence in the ratio of the number of old 
stars to the number of young stars. In contrast, 
no variation is seen in the [3.6] — [4.5] color, and 
these colors are consistent with old stars. Thus 
these two shorter wavelength bands are dominated 
by an older stellar population. Only a small vari- 
ation of 0.2 magnitudes is seen in the [5.8] — [8.0] 
color, and these colors are consistent with those of 
interstellar dust, thus these bands are dominated 
by interstellar matter. The [5.8] — [8.0] color be- 
comes redder in regions with stronger star forma- 
tion, as indicated by the Ha and mid-infrared lu- 
minosities. 

Our numerical model suggests that Arp 107 is 
the result of a hybrid collision: part small impact 
parameter/large angle of attack, as in colliding 
ring galaxies, and part prograde planar encounter, 
as in M51-type bridge/tail systems. Specifically, 
the model shows that the strongest compressions 
occur at different times along the primary arm. 
Assuming these compressions trigger star forma- 
tion, this differential delay can account for the ob- 
served variation of clump age with azimuth. The 
Spitzer observations have allowed us to detect this 
trend, and the model comparison allows us to un- 
derstand the driving mechanism. This example 
suggests bright prospects for learning about star 
formation mechanisms by comparing Spitzer data 
to models of interacting galaxies. 

We thank the Spitzer team for making this re- 
search possible. We also thank Bill Keel for pro- 
viding us with the velocity of the companion and 
the spectral type of the foreground star, and Mar- 
tin Cohen for providing stellar colors. We appreci- 
ate helpful suggestions on this project from Mark 
Giroux and the referee Suzanne Madden. This 
research was supported by NASA Spitzer grant 
1263924 and NSF grant AST-0097616. This re- 
search has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by 
the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration. 
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Fig. 1. — The Pick R band image of Arp 107. Note 
the strong primary arm/partial ring in the south, the 
tail to the northwest, the bridge, and the plume to the 
northeast of the companion. 



Fig. 2. — The continuum-subtracted Hq-|-[N II] im- 
age of Arp 107. Several star forming regions in the 
ring/arm are detected, along with the Seyfert nucleus. 
Some artifacts are present near the companion due to 
imperfect continuum subtraction. 
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Fig. 3. — Left to right, top to bottom: the Spitzer 
3.6, 4.5, 5.8, 8.0, and 24 urn maps of Arp 107. The 
field of view is the same as in Figures 1 and 2. 



Fig. 4. — The Spitzer 8.0 fim map of Arp 107, with 5" 
radius circles identifying the clumps discussed in the 
text. The clumps are numbered in order by increasing 
declination (see Table 1). The number identifying each 
clump is directly above the corresponding circle in this 
figure. 



Fig. 5.— The Spitzer IRAC [3.6] - [4.5] vs. [4.5] - 
[5.8] color-color plot for the brightest 8 /xm clumps in 
Arp 107 (black crosses). The clumps are labeled. The 
colors of MOIII stars (open dark blue square), from M. 
Cohen (2005, private communication), and the mean 
colors of the field stars of Whitney et al. (2004) (ma- 
genta triangle) are also shown. The colors of normal 
stars all lie within 0.5 magnitudes of 0, in this plot 
(M. Cohen 2005, private communication). The Hatz- 
iminaoglou et al. (2005) colors of quasars are also plot- 
ted (red circles), along with the colors of the nuclei 
of the two galaxies (#14 and #28). The predicted 
IRAC colors for interstellar dust (Li & Draine 2001) 
are also plotted (green diamonds), for ISRF strengths 
that vary from 0.3 — 10,000 x that in the solar neigh- 
borhood. As the radiation field increases, the colors 
become redder. Note that the predicted dust colors 
vary very little for this wide range in ISRF. The er- 
rorbars include both statistical uncertainties and an 
uncertainty in the colors due to varying the sky annuli 
(see text). 
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Fig. 6.— The Spitzer IRAC [4.5] - [5.8] vs. [5.8] - 
[8.0] color-color plot for the brightest 8 fim clumps in 
Arp 107 (black crosses). The clumps are labeled. The 
colors of MOIII stars (open dark blue square), from M. 
Cohen (2005, private communication), and the mean 
colors of the field stars of Whitney et al. (2004) (ma- 
genta triangle) are also shown. The colors of normal 
stars all lie within 0.5 magnitudes of 0, in this plot 
(M. Cohen 2005, private communication). The Hatz- 
iminaoglou et al. (2005) colors of quasars are also plot- 
ted (red circles), along with the colors of the nuclei 
of the two galaxies (#14 and #28). The predicted 
IRAC colors for interstellar dust (Li & Draine 2001) 
are also plotted (green diamonds), for ISRF strengths 
that vary from 0.3 — 10,000 x that in the solar neigh- 
borhood. As the ISRF increases, the IRAC colors of 
dust become redder. Note that the predicted dust col- 
ors vary very little for this wide range in ISRF. The 
errorbars include both statistical uncertainties and an 
uncertainty in the colors due to varying the sky annuli 
(see text). 
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Fig. 7.— The Spitzer IRAC [3.6] - [4.5] vs. [4.5] 
— [5.8] color-color plot for the fainter 8 fim clumps in 
Figure 4 (black crosses) . The clumps are labeled. The 
colors of MOIII stars (open dark blue square), from M. 
Cohen (2005, private communication), the mean colors 
of the field stars of Whitney et al. (2004) (magenta tri- 
angle), the predicted IRAC colors for interstellar dust 
(Li & Draine 2001) (green diamonds), and the Hatz- 
iminaoglou et al. (2005) quasars colors (red circles) are 
also plotted. The very faint sources #18 and #25 are 
not plotted since their uncertainties are very large and 
their colors are not well-constrained. 
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Fig. 8.— The Spitzer IRAC [4.5] - [5.8] vs. [5.8] 
— [8.0] color-color plot for the fainter 8 /im clumps 
in Figure 4 (black crosses). Some of the clumps are 
labeled. The colors of MOIII stars (open dark blue 
square), from M. Cohen (2005, private communica- 
tion), the mean colors of the field stars of Whitney et 
al. (2004) (magenta triangle), the predicted IRAC col- 
ors for interstellar dust (Li & Draine 2001) (green dia- 
monds), and the Hatziminaoglou et al. (2005) quasars 
colors (red circles) are also plotted. The very faint 
sources #18, #22, and #25 are not plotted since their 
uncertainties are very large and their colors are not 
well-constrained. 
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Fig. 9. — The Spitzer IRAC colors for the clumps in 
the ring/arm, as a function of position angle. Clumps 
#20 and #22 are not included in the second and third 
panels since the 5.8 /im uncertainties are very large. 
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Fig. 11. — For the clumps in the ring/arm, the Spitzer 
[8.0] — [24] color as a function of position angle. For 
clarity, clumps #15, 17, 20, and 22 are not plotted, 
since their errorbars are very large. 
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Fig. 10. — For the identified clumps, the Spitzer 
IRAC [5.8] - [8.0] color plotted against the IRAC- 
MIPS [8.0] - [24] color (black crosses). The clumps 
are labeled to the upper left of the corresponding data 
point. The colors of the Hatziminaoglou et al. (2005) 
quasars colors (red circles) are also plotted, along with 
the predicted colors for interstellar dust (Li & Draine 
2001) (green diamonds) . The green labels on the data 
points for dust correspond to the ISRF in terms of that 
in the solar neighborhood. The errorbars include both 
statistical uncertainties as well as an extra term due to 
background region selection, as discussed in the text. 
For clarity, clumps #15, 17, 20, and 22 are omitted 
because they have very large errorbars. 
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Fig. 12. — Top: For the identified clumps, the ratio 
of Ha luminosity to 3.6 /im luminosity, as a function 
of 8.0 /im magnitude. Some of the clumps are labeled. 
Bottom: the ratio of Hq luminosity to 8.0 /im lumi- 
nosity, as a function of 8.0 magnitude. The upper 
limits are 3a. 
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Fig. 13. — Top: the ratio of Ho luminosity to 3.6 
fim luminosity as a function of position angle, for the 
clumps in the ring/arm. Bottom: the ratio of Ha lu- 
minosity to 8.0 luminosity as a function of position 
angle, for the clumps in the ring/arm. The plotted er- 
rorbars only include statistical errors, not uncertain- 
ties in the absolute photometry. 



Fig. 14. — Four snapshots from a numerical hydro- 
dynamical model of the Arp 107 system. The top left 
panel shows a sky plane view of gas and star particles 
at a time near the present. The dashed line shows the 
relative orbit of the companion from the time of the 
first close approach. The top right panel is the same, 
but viewed from an orthogonal perspective. The lower 
panels show curves of test particles that were circular 
rings in the initial primary disk, at two earlier times, 
in the sky plane. In both cases the extension to the 
west is made up of particles in the bridge or captured 
onto the companion. Those two panels show the de- 
velopment of compression at the front of the primary 
arm, which is the likely driver of the observed star for- 
mation in this system. Length units are in kpc, with 
the adopted scaling. See text for details. 
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Table 1 

Spitzer Clumps in the Arp 107 Field 



ID 




R.A. 


(J2000) 


Dec. 




F3.6/<m 

(mJy) 


F4.5nTn 

(mJy) 


(mJy) 


(mJy) 


F24,»m 

(mJy) 


^Hcx + [N II] 

(10-1= 
erg cm-^) 


1 


10 


52 


15.8 


30 


2 


42.9 


0.123 


0.078 


0.045 


0.086 


0.461 


<1.0 


2 


10 


52 


14.0 


30 


2 


43.5 


0.079 


0.056 


0.057 


0.140 


<0.044 


<1.0 


3 


10 


52 


17.8 


30 


2 


51.3 


0.132 


0.103 


0.068 


0.159 


0.258 


<1.0 


4 


10 


52 


15.4 


30 


2 


56.7 


0.279 


0.176 


0.417 


1.367 


0.907 


1.7 


5 


10 


52 


14.7 


30 


2 


58.1 


0.237 


0.159 


0.375 


1.268 


0.979 


2.7 


6 


10 


52 


15.9 


30 


3 


0.2 


0.289 


0.177 


0.313 


0.957 


0.613 


<1.0 


7 


10 


52 


13.6 


30 


3 


2.6 


0.207 


0.128 


0.350 


1.001 


0.584 


2.6 


8 


10 


52 


17.0 


30 


3 


11.7 


0.155 


0.095 


0.119 


0.548 


0.227 


<1.0 


9 


10 


52 


14.1 


30 


3 


12.7 


1.685 


1.079 


0.732 


0.336 


<0.075 


<1.0 


10 


10 


52 


12.8 


30 


3 


17.3 


0.180 


0.127 


0.458 


1.478 


1.733 


12.8 


11 


10 


52 


16.3 


30 


3 


18.3 


0.057 


0.027 


0.027 


0.111 


0.114 


<1.0 


12 


10 


52 


18.6 


30 


3 


21.2 


0.306 


0.209 


0.162 


0.177 


0.546 


<1.0 


13 


10 


52 


14.0 


30 


3 


26.7 


0.105 


0.061 


0.088 


0.202 


1.824 


<1.0 


14 


10 


52 


15.0 


30 


3 


28.4 


2.793 


1.848 


2.003 


4.250 


36.603 


<1.0 


15 


10 


52 


17.2 


30 


3 


29.0 


0.146 


0.083 


0.174 


0.528 


0.584 


<1.0 


16 


10 


52 


12.6 


30 


3 


32.3 


0.042 


0.026 


0.059 


0.269 


0.486 


3.6 


17 


10 


52 


16.9 


30 


3 


37.7 


0.157 


0.088 


0.125 


0.406 


0.679 


<1.0 


18 


10 


52 


13.6 


30 


3 


39.0 


0.005 


<0.001 


< 0.005 


0.053 


0.258 


<1.0 


19 


10 


52 


10.9 


30 


3 


45.1 


0.038 


0.039 


0.067 


0.112 


0.100 


<1.0 


20 


10 


52 


16.2 


30 


3 


47.1 


0.071 


0.040 


<0.006 


0.064 


0.290 


<1.0 


21 


10 


52 


12.8 


30 


3 


49.1 


0.130 


0.081 


0.348 


1.298 


1.087 


5.8 


22 


10 


52 


15.7 


30 


3 


54.3 


0.060 


0.042 


0.015 


0.195 


0.446 


<1.0 


23 


10 


52 


9.9 


30 


3 


56.8 


0.147 


0.205 


0.256 


0.361 


1.219 


<1.0 


24 


10 


52 


15.8 


30 


4 


1.1 


0.054 


0.038 


0.061 


0.224 


0.443 


<1.0 


25 


10 


52 


11.3 


30 


4 


1.5 


0.133 


0.092 


< 0.005 


0.037 


<0.041 


<1.0 


26 


10 


52 


14.7 


30 


4 


5.3 


0.188 


0.116 


0.517 


1.913 


1.760 


8.7 


27 


10 


52 


17.6 


30 


4 


8.1 


0.116 


0.097 


0.093 


0.141 


0.159 


<1.0 


28 


10 


52 


18.5 


30 


4 


20.9 


5.457 


3.145 


2.326 


1.530 


0.546 


<1.0 


29 


10 


52 


12.6 


30 


4 


21.3 


0.072 


0.042 


0.039 


0.190 


0.122 


<1.0 



"^The statistical uncertainties at 3.6, 4.5, 5.8, 8.0, and 24 jum are <0.002, <0.001, <0.002, <0.005, and <0.013 mJy, 
respectively, determined from the rms in the sky annuH. The absolute calibration uncertainty is <10% (IRAC Data 
Handbook). 

^The statistical Ha+[N II] uncertainties are paS.S X 10~^^ erg s~^ cm~^. The absolute Ha+[N II] calibration 
uncertainty is Ri30%. 
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